Evolution of linear perturbations through a bouncing world model: 
Is the near Harrison-Zel'dovich spectrum possible via a bounce? 
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We present a detailed numerical study of the evolutions of cosmological linear perturbations 
through a simple bouncing world model based on two scalar fields. We properly identify the rela- 
tively growing and decaying solutions in expanding and collapsing phases. Using a decomposition 
based on the large-scale limit exact solution of curvature (adiabatic) perturbations with two indepen- 
dent modes, we assign the relatively growing/decaying one in an expanding phase as the C/d- mode. 
In the collapsing phase, the roles are reversed, and the C/d- mode is relatively decaying/growing. 
The analytic solution shows that, as long as the large scale and the adiabatic conditions are met, 
the C- and d-modes preserve their nature throughout the bounce. Here, by using a concrete non- 
singular bouncing world model based on two scalar fields, we numerically follow the evolutions 
of the correctly identified C- and d-modes which preserve their nature through the bounce, thus 
confirming our previous anticipation based on the analytic solution. Since we are currently in an 
expanding phase the observationally relevant one in the expanding phase is the relatively growing 
C-mode whose nature is preserved throughout the bounce. The spectrum of C-mode generated 
from quantum fluctuations in the collapsing phase has a quite blue spectrum compared with the 
near Harrison-Zel'dovich scale-invariant one. Thus, while the large-scale condition is satisfied and 
the adiabatic condition is met during the bounce, we conclude that it is not possible to obtain the 
near Harrison-Zel'dovich scale-invariant density spectrum through a bouncing world model as long 
as the seed fluctuations were generated from quantum fluctuations of the curvature perturbation in 
the collapsing phase. We also study the tensor-type perturbation. For the tensor-type perturbation, 
however, both C- and d-modes (of the tensor-type perturbations) in the collapsing phase survive as 
the relatively growing C-mode in the expanding phase. Due to its growing nature in the collapsing 
phase, the initial d-mode dominates the surviving C-mode spectrum after the bounce. 

PACS numbers: 98.80.-k, 98.80.Cq, 98.80.Jk 
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I. INTRODUCTION 

In order to explain the large-scale structures and 
motions of the galaxy distribution and the tempera- 
ture anisotropics of cosmic microwave background we 
need near Harrison-Zel'dovich scale-invariant spectrum 
for the initial density fluctuations which was suggested 
in the early 70's Currently, inflation (early accel- 
eration phase) scenario provides the only known mech- 
anism which can successfully generate seed fluctuations 
for such cosmological structures. Inflation provides near 
Harrison-Zel'dovich spectrum rather naturally by ampli- 
fying the ever present quantum fluctuations in the matter 
and space-time in microscopic scales during an early ac- 
celeration phase. 

Recently, several bouncing cosmological world models 
were suggested as alternatives to the inflation scenario. 
The pre-big bang scenario @ is based on the low-energy 
tree-level string theory effective action. However, pre- 
big bang scenario fails to explain the current large-scale 
structure because the generated spectrum is quite bluer 
than near Harrison-Zel'dovich spectrum, see [|| . More re- 
cently, ekpyrotic and cyclic scenarios [J] were suggested 



2 



based on the collision of two branes in a five dimensional 
bulk. According to the authors of these scenarios, near 
Harrison-Zel'dovich spectrum follows from the d-mode 
of curvature perturbation in the zero-shear gauge (y x ; 
terms will be explained shortly). Such a claim has been 
disputed by several authors [EE UMM- We Pointed out 
that although the d-mode is the dominating mode in the 
collapsing phase it decays away as the background model 
switches to expanding phase, thus irrelevant observation- 
ally 0, H, [H . Also the d-mode of ip x is known to exag- 
gerate the perturbations due to overly distorted hyper- 
surface (temporal gauge) condition compared with the 
curvature perturbation in the comoving gauge (tp v ), see 
[Icj |. In fact, the authors of cyclic scenario in [H claimed 
that in effective four dimensional space-time their world 
model goes through a singularity. In such a case not 
only the background world model, but also the pertur- 
bations inevitably become singular, thus we have no tool 
to handle the situation 5]. In order to handle the situa- 
tion theoretically in @, d, Q we have considered several 
smooth and nonsingular bouncing world model assuming 
both the background world model and linear perturba- 
tions remain valid. The present work is a continuation 
of our previous studies using a specific bouncing world 
model and concrete numerical study of perturbations un- 
der such a bounce. 

Here we summarize our previous reason, for not consid- 
ering the d-mode as seed fluctuations, used in 0,1,1. a 

scalar-type cosmological perturbation of a single compo- 
nent medium is described by a second-order differential 
equation. As long as the large-scale condition is met we 
have an analytic solution with two independent modes 
which we term C- and d-modes, see Sec. IIII Al We as- 
sign the relatively growing/decaying one in an expand- 
ing phase as the C /d-mode; in the collapsing phase, the 
roles are reversed, and the C/d-mode is relatively decay- 
ing/growing. In fact, for tp v the C-mode simply remains 
constant in time and the d-mode decays in an expanding 
phase, thus suggesting their names. Since such a large- 
scale condition is well met near the bounce we are con- 
sidering, we anticipate the nature of the C- and <i-modes 
to be preserved even after the bounce. Because the ob- 
servationally relevant one in the expanding phase is the 
relatively growing C-mode, even in the collapsing phase, 
we have to consider the seed perturbation in the C-mode 
generated from quantum fluctuations. This shows quite 
blue spectra both for the pre-big bang and the ekpy- 
rotic (or cyclic) scenarios. Although not successful as 
alternative to the inflation, these models have spurred 
interests in the evolution of perturbations through non- 
singular bounces in the early universe [Til . H3 , Il3| . In 
this work we will show the same conclusion, now using a 
numerical study of perturbation through a specific non- 
singular bouncing world model. 

In order to have a realistic bounce in practice, we have 
to consider an additional field or fluid component which 
drives the nonsingular bounce in additional to the field 
or fluid which gives rise to the seed fluctuations for the 



large-scale structures [9|, In such cases we have to 
consider perturbations of the additional field introduced 
to have bounce, in addition to the field for fluctuation 
seeds. Even in such cases, however, we expect our above 
conclusion based on the analytic solutions for the single 
component is relevant for the analysis; for our previous 
analytic argument, see Sec. V.C of [9[. Recently, however, 
authors of [12| have reported that by numerical analysis, 
for a specific collapse phase, the d-mode of ip v in the col- 
lapsing phase survives as the C-mode in the expanding 
phase, thus allowing a scale-invariant spectrum available 
for that specific collapse, see also @, [IH, The nu- 

merical results in [12| show the growing mode (d-mode) 
in a collapsing phase survives as the growing mode (C- 
mode) in an expanding phase. In this work we will show 
that this is a result based on mixed initial conditions 
of our precise C- and d-modes; the latter initial condi- 
tions produce evolutions where C- and d-mode natures 
are preserved even after the bounce. 

In this work, we would like to address and clarify 
whether switching of d-mode to the C-mode is possible 
as the perturbation goes through a nonsingular bounce. 
We will resolve such an issue by numerically following 
the evolution of diverse C- and d-mode initial conditions 
through a simple toy bounce model based on two scalar 
fields. In order to have the C- and d-mode decomposi- 
tion, the the curvature (adiabatic) perturbations should 
evolve independently (we will call it the adiabatic condi- 
tion) from the isocurvature perturbations, and the large- 
scale condition should be met. Our results show that, 
as long as the large-scale condition and adiabatic con- 
dition are met, we can find our C- and d-mode initial 
conditions where the nature of C- and d-mode are pre- 
served through the bounce; as some background variables 
diverge near the bounce we will encounter some sharp 
breakdowns of the above conditions caused by the di- 
vergence of the background variables; we will argue that 
such sharp breakdowns do not have physical impacts on 
interpreting the C- and d-mode natures. 

In order to get these results which are supported by 
previous analytic methods in Q , it is important to find 
the precise C- and d-mode initial conditions which allow 
the proper evolutions of these two modes. Since our real- 
istic bounce model is based on two scalar fields, the naive 
analytic form C- and d-mode initial conditions based on 
analytic solutions of a single component medium are not 
the precise ones. We will show that the previous dif- 
ferent reports were based on taking these naive initial 
conditions which in fact correspond to mixed initial con- 
ditions of our (precise) C- and d-modes. Under such 
mixed initial conditions the adiabatic condition breaks 
down near the bounce, thus one cannot trace the re- 
sulting C-mode in the expanding phase to the d-mode 
in the collapsing phase. We will show that under such 
mixed initial conditions the isocurvature perturbations 
are excited near the bounce. From our study we will 
conclude that, as long as the large-scale conditions and 
the adiabatic conditions are met, it is not possible to have 
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near Harrison-Zel'dovich spectrum from the C-mode gen- 
erated from quantum fluctuations during the collapsing 
phase before bounce. We emphasize the conditions we 
use to get such a strong conclusion: (i) Einstein's gravity 
and linear perturbation theory work throughout the evo- 
lution, (ii) quantum fluctuations of a minimally coupled 
scalar field in a collapsing phase provide seed curvature 
(adiabatic) mode fluctuations for the large-scale struc- 
ture, (iii) the large-scale conditions are met during the 
bounce, (iv) the adiabatic conditions are met during the 
bounce. We will show that the latter two conditions are 
often violated subtly near the bounce without affecting 
curvature perturbation natures in the large-scale limit. 

In the case of the tensor-type perturbation (gravita- 
tional waves) we will present the C-mode which remains 
constant throughout the bounce. However, the d-mode 
in the collapsing phase generates an additional C-mode 
after the bounce. Since the d-mode grows very rapidly 
during the collapsing phase, the C-mode generated by the 
initial d-mode is likely to be dominant over the original 
C-mode initially generated during the collapsing phase. 
In such a case, it should be the d-mode initial condition 
in the collapsing phase which gives rise to the dominat- 
ing C-mode tensor-type perturbation in an expanding 
phase after the bounce. We will show that for a collaps- 
ing phase with matter-dominated-like equation of state 
we could have near Harrison-Zel'dovich spectrum for the 
tensor-type perturbation. Its amplitude, however, should 
depend on the duration of the collapsing phase after its 
generation from quantum fluctuations. 

In Sec. |H]we introduce our nonsingular bouncing world 
model based on an ordinary minimally coupled scalar 
field together with a ghost scalar field without potential; 
the latter field causes a nonsingular bounce. In Sec. lIIII wc 
present basic perturbation equations, large-scale analytic 
solutions and exact solutions under special situations in 
three different temporal gauge conditions. We include 
the tensor- type perturbations. In Sec. IIVI we present 
our numerical integrations based on our precise initial 
conditions which give correct behaviors of the C- and d- 
modes in both collapsing and expanding phases. We will 
check the conditions required to have such a decompo- 
sition. We compare our results with the ones based on 
the naive initial conditions which are motivated by the 
analytic initial conditions in a single component situa- 
tion. We will show that these latter initial conditions 
in fact correspond to the mixed initial conditions, and 
show that as the evolution proceeds certain conditions 
needed for the mode decomposition are violated near the 
bounce. We also analyze the tensor-type perturbations. 
Since perturbation equations and variables often show 
singular behaviors near the bounce, in this section we 
carefully examine and explain the evolutions and the re- 
quired conditions. We also explain the method of finding 
our precise C- and d-mode initial conditions. In Sec. [V] 
we consider the quantum fluctuations as the initial con- 
ditions and present the generated power spectra for the 
scalar- and tensor- type perturbations. We also present 



implications of our results. Section IVII presents discus- 
sion. 

In this work, we assume Einstein's gravity and spa- 
tially homogeneous-isotropic background world model 
with linear perturbations valid throughout the bounce. 
We set 8ttG = 1 = c and H = 1. 



II. BACKGROUND 

In order to have a collapsing phase followed by an ex- 
panding phase which is smoothly connected by a non- 
singular bounce we consider a model based on two scalar 
fields introduced by Allen and Wands in [l^] • One field </> 
is an ordinary minimally coupled scalar field with a pos- 
itive kinetic energy and a self-interaction potential. The 
other field a has a negative kinetic energy (a ghost field) 
without a potential. The ghost field could violates null- 
energy condition [l5|, [i + p < 0, which is necessary to 
have a smooth and nonsingular bounce. The ghost scalar 
field dominates near the singularity and causes a non- 
singular bounce. The action and the energy-momentum 
tensor are 



S = 



T„ 



d x. 



f 



dab 



(1) 



(2) 



where tildes indicate covariant variables; the Latin in- 
dices indicate the space-time. 

As the background world model we consider a spatially 
flat, homogeneous and isotropic, Friedmann model. To 
the background order the Robertson- Walker metric can 
be written as 



ds 2 = 



a 2 (-dry 2 



{3) dx a d X P 



9 a0 



(3) 



where a(rf) is the cosmic scale factor, rj is a conformal 
time, and g a i could become S a p in our flat background; 
the Greek indices indicate the space. The effective fluid 
quantities are 



f 



V, p 



1 



1 



V, (4) 



where p, and p are the energy density and pressure, re- 
spectively; a dot indicates a time derivative based on the 
time t with dt = adr]. The basic equations for background 



B* = i ( '-<P 2 - -a 2 
3 V2 2 



V 



H 



-5< 

3#0- 
3Ha 



o, 



0. 



(5) 

(6) 

(7) 
(8) 
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where H = a/a. Equation ([5]) and ([6]) are the Friedmann 
equations, and Eqs. © and §8§ are Klein-Gordon equa- 
tions for each of the scalar fields. We take the potential of 
the ordinary scalar field to be a simple exponential form 



V = V cxp(-X(j)). 



(9) 



We solve Eqs. ©-(O numerically to follow the back- 
ground evolution. We can integrate Eq. ([5]) to obtain 



a oc 



Thus, the energy density and pressure of the ghost field 
dominate near the bounce. Figure [T] shows the evolution 
of scale factor and its derivative near the bounce. Fig- 
ure [2] shows that the energy density of ghost scalar field 
dominates only near the bounce. Away from the bounce, 
where the ghost scalar field is negligible, assuming the 
exponential potential in Eq. ^ we have the power-law 
expansion or collapse pjj 

o o 

a oc t q oc r\ 1 -i , q 



3(1 + ^) A 2: 



V 



V 3 ( 1+ " 



= J3(l + w<,)H, (11) 



where 

b 2 - V. Thus, we have A = ^3(1 



Pw/fJ>W witn H4>) = \v + V and P(<P) = 

Wrf,). The value of 



A determines the evolution of scale factor far away from 
the bounce. For A = 2 we have w$ = 1/3 and the scale 
factor evolves like the radiation dominate era, and for 
A = \/3 we have w$ = which is a matter dominate era. 
In this work we will take A = y/3, thus w<p — 0. 

In Figs. Q] and [H we show the background evolution 
through the bounce. In order to help numerical repro- 
duction, we present the precise numerical value of our 
initial conditions in Table |TJ A prime indicates a time 
derivative based on the conformal time rj. 



TABLE I: The initial conditions we used for the background 
evolution. We set a m i n = 1, Vb = 0.01, and A = a/3- 



variable 



initial value 



r, -5.0E+02 

1.05727052142872E+01 

<j>' -6.6836661819519E-03 

a -1.8136774745271E+00 

a' 7.0552762997745E-07 

a 4.4771380493724E+02 

a! -1.7276454036056E+00 



III. PERTURBATIONS 

We introduce completely general perturbed order vari- 
ables in the Robertson- Walker background. The metric 



becomes 

5oo = -a 2 (1 + 2a) 



9o* = -a 2 (j3, a + , 

(12) 

where an index indicates 77; a vertical bar | indi- 



es) 

cates a covariant derivative based on g a L , and 2C ( 



(10) <#> + c 



("1/3) - 



|/3 T ^/8|a' perturbed order variables a, (3, 7 and 
if indicate the scalar- type perturbations. The transverse 

(v) (v) 

B a and C a indicate the vector- type perturbations. 
The transverse-tracefree C^l indicates the tensor-type 
perturbation. Indices (v) and (t) indicate the vector- 
and tensor-type perturbations, respectively. In our back- 
ground world model the three types of perturbations 
evolve independently to the linear order. We follow con- 
ventions in (l7l [Til . Il9j . The energy-momentum tensor 
can be decomposed into the effective fluid quantities as 

T ° = -(/x + , T Q ° = ( Ji + p)(-v, a + «W ) , 
T$ = (p + 6p)S% + W$, (13) 

where tracefree ITS indicates anisotropic stress. For the 
scalar fields we set 



a + 5a. 



(14) 



The two scalar fields contribute to the effective fluid 
quantities as 



5/j, - 
Sp = 4>Scf> - 



{n+p)v = - 
a 

11% =0. 



- 6S& 

&S& - 
1 



<r 2 )a - 



&Sa), (ji + p)vW=0, 



a 2 )a - 



(15) 



When we solve the Einstein equations, we have free- 
dom to choose the spatial and temporal gauge conditions. 
Practically, it is important to take a gauge which suits 
the problem, but generally we do not know the suitable 
gauge condition, a priori. The tensor-type perturbation 
is gauge-invariant and the vector-type perturbation can 
be written in a uniquely gauge-invariant form. By using 
a spatially gauge-invariant combination 



X = a(/3 + a A f), 



(16) 



instead of (3 and 7 individually, all scalar-type pertur- 
bation variables are spatially gauge-invariant. For the 
scalar-type perturbation, a proposal made in (TtI [Lsj 
is that we write the set of equation without fixing the 
temporal gauge (hypersurface) condition and arrange the 
equations so that we can implement easily various funda- 
mental temporal gauge conditions. The temporal gauge 
conditions should be chosen to achieve either mathemat- 
ical simplicity or physical interpretation. 
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FIG. 1: Evolution of the scale factor and its time derivative. In all figures in this work we take A = \fi. This leads to 
matter-dominated-era-like behaviors (a oc rf oc 

t 2/ 3) 

away from the bounce. Panel (a) shows evolution of the scale factor. The 
bounce is nonsingular, and we set a min = 1 at the bounce. Panel (b) shows time derivative of the scale factor. At the bounce, 
the derivative of the scale factor becomes zero. 



-20 



-40 



-60 - 



-80 - 



100 



_J I L_ 



-2 



(a) 



sign(cO -log(cv/cv m J 



o 0. 
is 



o 
o 



0.4 



u 

CD 

c 



(b) 



~i 1 r 



Im ( „)/m, 



(0)1 



-2 




sign(cO -log(a/a m J 



FIG. 2: Evolutions of the equation of state and the energy density ratio of the two fields. Panel (a) shows evolution of the 
equation of state w = The equation of state w is zero at early and late times which behaves like the matter dominated 

eras. Panel (b) shows the energy density ratio of the two scalar fields: fj,^ = + V and /U( CT ) = — \a 1 . The energy density 
of ordinary scalar field (j> dominates except near the bounce where the ghost field a becomes significant. The bounce is caused 
by the presence of a. The horizontal axis indicates the time axis using the scale factor; time increases from left to right, and 
— (+) sign denotes the collapsing (expanding) phase. 
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The complete set of equations describing the scalar- 
type perturbation without fixing the temporal gauge con- 
dition, i.e., in a gauge-ready form, can be found in Eqs. 



(10)-(15) of [IJ. Using Eqs. gj and £5) for the fluid 
quantities we have 



k = 3 (-ip + Ha) - ^?x, 



A 1 

—rip + Hn = 



— aSa — 



a* ■ 2 

K + -^x = g (<W _ &6<7 ) > 

X + Hx~a-(p = 0, 

k + 2Hn = - + 3H^j a + |l 
A 



- cr 2 )a + T^ <50 



(/>50 — (7(5(7 



6 2 -a 2 )a 



U <5</>} 



<5</> + 3if(50 + V^cty -54> = cj)(K + a+ 3Ha) + 20a, 

a 2 

5a + 3H5a -5a = &(k + a + 2>Ha) + 2ira. 



(17) 

(18) 

(19) 
(20) 

(21) 
(22) 
(23) 



The scalar fields do not support the vector-type pertur- 
bations. The tensor- type perturbation is described by 

m 



a 



(*) 



3HC 



(*) 



0. 



(24) 



The presence of scalar fields affects the tensor-type per- 
turbation only through the background evolution. 

The gauge transformation property can be found in 
[I?], HH • Under the gauge transformation x a = x a + £ a 
with £' = a£° (0 = 77) and £ Q = ^ a + we have 

a = a-i\ <p = (p-H^, 
A 



Sfi = 5p — /i£ , 5p = 5p — pd, , v = v £ 

a 



</>f, 6a = 5a- <rf 



(25) 



Thus, all variables in our Eqs. (ri"7|) - (l2"3"]) are spatially 
gauge- invariant. It will be convenient to introduce sev- 
eral gauge invariant combinations 



<p v =ip- aHv, p x = ip- Hx, 
5a x = 5a — ax, 5(j)s a = 5(f> — -, 



5(f> - cf>x, 



-6a. 



(26) 



ip v is the same as ip in the comoving gauge which sets 
v = 0, etc. It also allows algebraic manipulations like 



<P64, = <f- 



H 



H 



(27) 



etc, which is convenient in practice. 

A. Large-scale solutions with C- and d-modes 

Here, we introduce the large-scale asymptotic solutions 
with the C- and d-mode decomposition. The evolution 
of a single component medium is described by a second- 
order differential equation which has two independent 
solutions. In the large-scale asymptotic limit, we de- 
compose these two solutions into the C-mode and the 
d-mode. Since we are considering two fields, we have to 
check whether and when such a decomposition is possi- 
ble. 

Using Eqs. (fT7 ]) -([2"0" ]) we have 



H 2 



(28) 



Using Eqs. jn||, (JTSJ), ((2DJ), ((5SJ), and ((5SJ) we can show 



k 2 H 



2H Tr • „ [54> 5a\ 

V^a 2 -f - — , (29) 



(m + p) 



Our notation of the gauge-invariant variables explicitly 
shows the variable and the gauge condition: for example, 



where k is the wave number in Fourier space with A = 
— k 2 . By combining Eqs. (|2"8"|) and (|2"9"|) we have 
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H 2 
(fi+p)a 3 



H 2 



(J-+P 
H 



H 2 ( a 



_a(n +p) \H 



k 2 

~o<Pv 



k 2 

*** 



H 



v" - I — - k 2 } v 



0?sj[i + p 

2H 2 

{fi+p)a 3 {(»+p)H 



(i/«) 



-k 2 



Sa 
a 



1 

M + P 



V^a 2 



Sa 
a 



(30) 



(31) 



where 



Zfy, 



zH 



<Px'> 



H 2 ' 



(32) 



We note that, despite the notation, z 2 defined in Eq. |32|) 
could have the negative value depending on the sign of 



(p +p). Indeed, in order to have a smooth nonsingular 
bounce it is necessary to change sign of (fx + p) near the 
bounce, see Fig. [3] 

In order to make Eqs. (j3"0|) and (f3"Tj) complete we need 
equations for (5<f>/cj) - 8a/a). From Eqs. US]), (HI), fl9]). 
and (|2"31 we can derive 



/' 


+ P 


a 3 1; 


(> 2 <7 2 


f8<t> Sa 




f> 2 (7 2 


M 


+ P 


\~J ^ 


l« 


+ P 


a 3 <; 


9a 2 


( S(j) Sa 


a 3 <j 


!> 2 a 2 




+ P 





-3H + 4 



V+P 



„2 



S(p Sa 



2 V„ 



-3i? - 



fc 2 

^,2 



£(7 
(7 



6 cr J H 



2 V^fc 2 

T2^X' 



(33) 
(34) 



where we used Eq. (|29|) to derive Eq. (f34|) . Equations 
(|3T))) , (|3"Tj) , (f3"3")) , and (fM)) provide a complete set of equa- 
tions in terms of the curvature (adiabatic) perturbation 
variable ip v or tp x , and the isocurvature perturbation vari- 
able (S(f>/4> — Safer). We introduce a notation for the 
isocurvature perturbation variable 



Sa 
a 



! >Sa 



(35) 



The above equations are known in the literature in a 
more general context with multiple minimally coupled 
scalar fields, see [2(J. Here we showed that the general 
equations in [2(| are satisfied when one of the scalar field 
has negative kinetic term. The gauge-invariant variables 
tp v and tp x are the spatial curvature perturbation vari- 
able if in the comoving gauge (v = 0) and the zero-shear 
gauge (x = 0), respectively. Thus, both of these two vari- 
ables represent the curvature (adiabatic) perturbation. 
The second-order differential equation for (Sa/a — 8(f) / '</») 
represents the isocurvature perturbation, see [2(j. In 
this work, we are mainly interested in the evolution of 
the curvature-type perturbation (we call it the curvature 
mode) where isocurvature perturbation is suppressed; we 
will show that even in such a case, the isocurvature per- 



turbation does not necessarily vanish, and needs to be 
monitored. 

The large-scale conditions are satisfied if k 2 terms are 
negligible in Eqs. (|3T))) and (|3"Tj) , i.e., 



(1/2)" 



k 2 « — , k 2 « 

z (1/z) 



(36) 



for f v and f x , respectively. The adiabatic conditions are 
satisfied if terms in the RHSs of Eqs. ([30]) and (f3i"j) are 
negligible compared with terms in the LHSs involving v 
and u, respectively. That is, our adiabatic conditions are 
satisfied if 



2H 2 



a{n + p) - k 2 ) ip v 
a 2 V^a 2 (f - £ ) 



< 1, 



(m + p) 



fx 



< 1, 



(37) 



for ip v and f x , respectively. In order to monitor these 
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conditions it is convenient to have 



z 
z 



(VfT 



-5o 2 V - \ + 2^F 2 - ^V^' 



+ - 



,' 2 - a' 2 



a d> — a 



aV + 3^- + — 

a z a' 



a 2 V.^' Z + 2aa'V, 4> 0' - a\V,^f 

A' 2 



+2— VV^' + 3a 4 (V^) 2 - 



2 a ,2 



(38) 



Later we will show that our large-scale conditions and 
adiabatic conditions are often sharply broken near the 
bounce. Such sharp divergences occur because some 



background variables like fi + p, z, H, and <j> often vanish 
or diverge near the bounce. Later we will argue that such 
a very brief timescale breaking of our conditions do not 
have physical impact on interpreting the evolution of cur- 
vature (adiabatic) perturbations in the large-scale limit. 
In Fig. [3] we present evolutions of z" /z, (l/z)"/(l/z), a', 
(j)', and (/i + p), for later use. 

If both the large scale condition and adiabatic condi- 
tion are satisfied, even in our two component situation, 
we have the following general solutions 



p„(M) = C(k)-2d{k)k 2 / -L, 



MM) = -C(k) a 



2 , V dr, 
HP* H 



(39) 



z 2 dr] + d{k) — . (40) 



Notice that the d-mode of ip v is (fc?7) 2 -higher order in the 
large-scale expansion compared with the d-mode of ip x . 
To the next order in the large-scale expansion we have 

M 



i Pv (k,T l )=C(k)il + k 



'§1* 



2 , Pd V 



2d(k)k 2 J V ^, 



V x {hn) = \c{k)^ j\ 2 d V + ^d(k)il + k 2 J" i (f* 3 *l) dr, - J\ 2 d v J" J | 



(41) 
(42) 



We note that the integrations in Eqs. (|39"|) - (|42"]) do not 
consider contributions from the lower bounds of integra- 
tions; this, in fact, becomes a subtle point in the case of 
the bounce. 

In an expanding universe, C-mode is relatively growing 
and d-mode is decaying. In a collapsing phase, however, 
the d-mode grows rapidly and the C-mode relatively de- 
cays. Actually, the C-mode usually remains constant in 
either phase. Since we are considering two fields, in order 
to use the concept of C- and d-mode decomposition, we 
have to check whether both (i) the adiabatic conditions 
and (ii) the large-scale conditions are satisfied through- 
out the evolution. 



B. Equations in three different gauges 

We will investigate the behaviors of C- and d-modes 
by following their evolutions from the collapsing phase 
to the expanding phase through our simple non-singular 
bounce model. In order to obtain numerical results, it is 
practically important to take a suitable temporal gauge 
condition. If we have the behavior of any one perturba- 
tion variable in any one gauge condition the rest of the 
variables even in the other gauge conditions can be alge- 



braically derived using our equations in the gauge-ready 
form in Eqs. (|17[1 - (|23|> and the gauge-invariant combina- 
tions; the latters are made from the gauge-transformation 
properties in Eq. (f2"5|) , as in Eq. (|26[) . We have cho- 
sen two gauge-invariant combinations ip v and ip x because 
of their conserved properties in the large-scale limit and 
popularity in the literature. We could not directly solve 
ip v , however, because its equations show singular behav- 
ior near the bounce. We prefer a gauge condition that 
has no singular behavior through the bounce. We find 
several suitable gauge conditions and numerically solve 
them independently The solution of any gauge-invariant 
variable should be the same independently of the gauge 
conditions. This naturally provides a way to check the 
numerical accuracy. For instance, (p v and (p x can be ob- 
tained from the solutions known in the uniform-er gauge 
(5a = 0) and in the zero-shear gauge (x = 0) as 

a 1 4>'84>$ a a' 4>'S(f) x — a'5a x 



fx = fs<7 - -^Xs<y- (43) 

In the following we consider three different temporal 
gauge conditions. In each of these temporal gauge con- 
ditions, the temporal gauge degree of freedom is fixed 
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FIG. 3: Evolutions of z'/z, (1/jz)" /(1/z), i-i + p near the bounce [Panel (a)], and <f> , a' [Panel (c)]. Panel (b) magnifies Panel 
(a) in vertical axis. In Panel (a), +p) is plotted with 500 times magnification. Panel (b) shows that z" j z = and fi + p — 
occur at different epochs; z" / z — diverges at fi + p = and at the bounce. These combinations appear in Eqs. (|30[1 - (|34| ) which 
we use to distinguish the adiabatic (curvature) and isocurvature modes, and the large and small scales. As these variables 
(except for a') diverge or vanish some of our large-scale conditions or adiabatic conditions often become singular for very brief 
periods, see Figs. IHIfJl and ITOHTSI 



completely, thus each of the variables in these gauge con- 
ditions has unique gauge- invariant counterparts, see [l8l |. 
Thus, we can regard all the variables as gauge-invariant 
ones. 



1. Uniform- curvature gauge 



The uniform-curvature gauge takes ip = as the tem- 
poral gauge condition. This gauge condition is particu- 
larly suitable to handle quantum fluctuations generated 
during collapsing phase which will provide the initial seed 
fluctuations for the large-scale structures. This provides 
our initial conditions for numerical integration. From 
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Eqs. ([T7 )) -([l9 l) .((25 ]) . and (J23J) we can derive 



b% + 2-5^ + kHct> v 



l Srt> ' 



a 2 „2, 



a' a' 



k use/, 



a' 



6a" + 2-Sa' + k 2 6a v - a 2 ^a' 2 V6a v + a 2 ^a' ( V<f>' + —V g 
a a' a' \ a 



a 

cl' „ , , o „ o a" /2, , ,- t a" //,,,#. O' 



a' \ a' v ' 



fc U5o 

a 



a 2 ,i 



-aV—^a usa + a—a V —<p + — usa 
a' a \ a * 1 



V, 



0, 



0,(44) 



(45) 



r 



where 



U$ 



aSa v , use/, = ai 



(46) 



Away from the bounce, the effect of a is negligible, and 
we have the power law regime. From Eq. (jlip we have 
4>' oc a' I a. Thus, we have 



k 2 - U S<j> = 0. (48) 

For a general A we have the following exact solutions 

= Vk\v\[ci(k)H^(k\ V \) + c 2 (k)Hi 2 \k\ v \)l 
3(w- 1) _ 3(7-1 



2(3iw + l) 2(g-l) : 



(49) 



where (i) = <50 or Sa; A is related to q and u> as in Eq. 
(fTTjl . If we take A = V3, thus u> = and q = 2/3, we 
have 



These will provide the initial seed fluctuations. Identi- 
fication of these classical solutions with the initial fluc- 
tuations generated from the quantum fluctuations will 
be considered in Sec. [Vj The equations in this gauge, 
however, show the singular behavior as the background 
passes through the bounce. Thus, in order to follow the 
evolution numerically we need to use other gauge condi- 
tions which are well-defined throughout the bounce. 



2. Uniform-a gauge 



'(') 



+i 



cos{\krj\) 



sm{\kr]\) 
sin(|£;?7|) 



cos(|fc?7|) 



[ci(fc)+ca(fc)] 



M 



[ci(fc)-c 2 (fc)] \. (50) 



The uniform-cr gauge takes 8a = as the temporal 
gauge condition. Authors of [HJ used this gauge con- 
dition for their numerical calculation. From Eqs. (JT7J) , 
(fT9| . (|20| . (|22l) . and ([23]) we have 



+ 2—5<f>' Sc + a 2 V M 8(i> Sa + k 2 5<p s<J + 2a 2 V tj) a 5cJ = 0, 



a 2 



a L - 3<pL h — XSa = o, 



l V :4> 54>sa) - a 2 Va Sa a' Sa . = 0, 



x'sa H XSa ~ a(a Sa + <ps a ) = 0. 



(51) 
(52) 
(53) 
(54) 



Thus, in this gauge we solve sixth-order differential equations. Although, this set of sixth-order differential equations 
can be easily reduced to a fourth-order differential equation, in order to follow numerically the evolution of perturba- 
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tions through the bounce without singular behavior these sixth-order differential equations are preferred. In order to 
provide initial conditions for ag^ and xsa we use the following constraint equations 



asa 



XSo 



a 
a/ 



3 1 

2fc" 2 



1 a 
Wa 1 



<j>'6<t>' Sa - {4> - a' )asa + a 2 V^5^>s a 



(55) 
(56) 



which follow from Eqs. (fl~7[) and (fl"9]) and Eqs. (fl~8|) and (fT9|) . respectively. Since the initial conditions are provided 
in the uniform-curvature gauge, in order to construct initial conditions for the variables in the uniform-cr gauge it is 
convenient to have 



:8a u 



J 8a 



I , 2 T7 S(T <P 



<p 



a' Sa u 



a' 



a' 8a' 



(57) 



a a a a a 

From Eq. (|39p we can read the initial conditions for the C- and d-modes of ip v . In order to have the C- and d-mode 
decomposition, it is necessary to satisfy the large-scale and the adiabatic conditions. From Eq. ([30)) the adiabatic 
condition implies 8a / a — 6<p/4>. From this, we construct C- and d-modcs of ip v and the related variables in the 
uniform-cr gauge. Thus, the initial conditions for the C-mode and d-mode are, respectively 

C-mode : 5<j>g a = 0, 8<j>' Scr = 0, ip s<T = C, ip' Scr = 0, 
d-mode : 8<j>& a = 0, 5<f>' Sa . = 0, ipsa = d, (p' Sc 



= - 3 -d. 



(58) 



r 



It is important to notice, however, that these are ini- 
tial conditions for the C- and d-modes in an ideal situa- 
tion when the large-scale and the adiabatic conditions are 
satisfied exactly. Since both of these conditions are not 
satisfied exactly in the two-component medium, we have 
to find precise C- and d-mode initial conditions which 
have correct behaviors in the initial phase and preferably 
throughout the bounce. 

In Sec. IIVI we will present evolution of perturbations 
under such precise initial conditions (we call these 'our 
initial C- and d-mode' initial conditions) which preserve 
the C- and d-mode natures throughout the bounce; we 
will show that, under this evolution, the isocurvature per- 
turbation is not excited near the bounce. We will also 
present the perturbation behaviors under the naive initial 
conditions based on the analytic solutions in an ideal sit- 



uation in Eq. (|5"8"j) (we call these 'analytic C- and d-mode' 
initial conditions) which fail to produce the desired re- 
sults. In short, the initial conditions in Eq. I|58p. in fact, 
correspond to mixtures of the C- and d-modes in the ini- 
tial phase, and lead to loss of these mode identities as 
the evolution proceeds through the bounce; under this 
evolution, the isocurvature perturbation is also excited 
near the bounce, see Sec. IIVI 



3. Zero-shear gauge 

The zero-shear gauge takes x = as the temporal 
gauge condition. From Eqs. ([T7j). ([15]). [[2D]). ([22]), and 
(|2"3")l we have 



v'x 



a 



a'8a x ) 



Xh 
2 f j, 



= 40' 



5a'' + 2-Sa' + k 2 Sa x = 4a' 

x a x 



b x - a'5a x ) . 

In order to construct the initial conditions in this gauge we need the following relations 



a 



1 



^Vx + \{<t>' 5( t>x ~ a ' Sa x) 



2a V^<p x , 



—XSa, 5a x = -—Xs<r, 5(f>' = 8(j)' Sa + 3-^0'xfc + aV^xsa Xs<t> 



8a' 



S^a'xsa 



-XScr, Vx = fSa 



(59) 
(60) 
(61) 



(62) 
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where 



XSa = — < fSa 



1 

2fc2 

2 



* - CT )<PSo 



a 2 V 6 + -A'V 



V8a 



a 

2a' 



(63) 



Relations in Eq. ([63]) follow from Eqs. l[5l | -([56 j) . Equa- 
tions in this gauge also behave well throughout the 
bounce. Numerical results in this gauge condition should 
coincide with the ones from the uniform-cr gauge for 
any gauge-invariant combination which provides a self- 
consistent numerical check. 



C. Tensor-type perturbations 

The tensor-type perturbations correspond to the grav- 
itational waves. It is described by Eq. ([24]) . The scalar 
fields do not directly contribute to the tensor-type per- 
turbations, i.e., n a/ 3 = 0. Equation ([2"4"[) becomes 



3Mt) 



a0 



k 2 



a 



= o, 



(64) 



Eq. (|64|) has an exact solution 



Cl 



afl (k)HV>(k\T,\) 



+c 2afj (k)H%(k\r 1 \) 

I3q-l 
v„ = -- 




2 1-q 

cos(|fcri|) 4 
sm(\ki]\) 
ciapik) - c 2a p(k) 



sin(|fc7j|) 
IN 

cos(|fcri|) 



(68) 



\krj\ 



(69) 



where v t = aC^i. Notice the similarity of this equation 
compared with Eq. (|30[) for ip v . In the large-scale limit 
we have 



C$(*> 77) = C a0 (k) + d aP {k) I ^dr,, (65) 



" 1 



where C a p and d a $ are integration constants. To the 
next order in the large-scale expansion we have 



c$(M) = c^(fc)|i + fc 2 



71 / 1 r i x 

<r I —d-q ) dr/ 



drj a dr] 



"V 2 

+ d a p(k){ I -^dr\ 



V / 1 
a 



■^drj ) [ I -Ur? ) dr, 



* t o r> i 

(/ / — dr; ) dr? 



r*7 



The above solutions can be compared with Eqs. (|39|) . (|4T 
for ip v . 
For 



■(66) 



a rf 



g(2g-l) 



(67) 



IV. NUMERICAL RESULTS 
A. Scalar-type perturbations 

1. Our C- and d-modes 

In Figs. Ufa), Hfb), El^a), and [6] (b), we present evolu- 
tions of 'our C- and d-modes'. These are based on our 
C- and d-mode initial conditions presented in Table [Til 
These figures apparently show that evolutions away from 
the bounce coincide with the known behaviors of the C- 
and d-modes based on the large-scale analytic solutions 
in Eqs. ([39)) and (|40f . These analytic solutions are valid 
when both the large-scale conditions and the adiabatic 
conditions are satisfied. 

In Fig. [5] we examine the large-scale conditions using 
Eq. (j36p . Figure 02a) apparently shows that the large- 
scale condition for (p v is sharply broken four times near 
the bounce at z"/z — 0, see Fig. EH Despite these sharp 
divergences in our large-scale condition, we argue that 
these do not have physical impact on interpreting the C- 
and d-mode natures of ip v in Figs. Ufa) and [Ha). This is 
because our large-scale condition is sharply broken as the 
z" I z term in Eq. (|30[) crosses zero, thus becomes smaller 
than /c 2 -term, for very brief time intervals, see Fig. [3] 
Figure [SJb) shows that the large-scale condition is well 
satisfied for tp x . The large-scale conditions involve only 
the background evolution and the wavenumber. Thus, 
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the same large-scale condition applies to both C- and 
(i-modes independently of the initial conditions for per- 
turbations. 

In Figs.[Uc),[4jd),[6]Jc), an dEtd) we examine adiabatic 
conditions using Eq. ([37]) . All these figures show that the 
adiabatic conditions are sharply broken near the bounce. 
These occur at z" / z — 0, fi + p — 0, or ip x — 0. Despite 
these sharp divergences in our adiabatic conditions, we 
argue that these do not have physical impact on inter- 
preting the C- and d-mode natures of ip v in Figs. Ufa) and 
Ufa). We have two reasons. First, ip v in Figs. BJa) and 
[6f a) shows that the C- and d-mode natures are preserved 
away from the bounce. Second, later in Figs. [T"2l and [TBI 
we will show that for our C- and d-modes the isocurva- 
ture perturbation is not excited near the bounce. 

Compared with the above sharp spikes in the adiabatic 
conditions, Fig. [(Jfc) shows that the adiabatic condition 
for tp v is severely broken near the bounce region. Remem- 
ber that the d-mode of ip v in Eq. (f3T))) is higher order in 
the large-scale expansion compared with the same mode 
of ip x in Eq. (|40|) ; this happens because the leading or- 
der large-scale solution of the <i-mode of ip v has canceled 
out, see Eqs. (|28p and (|4"2")) . Thus, as the d-mode of ip v is 
(fc/ai/) 2 -factor smaller than the d-mode of tp x we could 
have such a breakdown in Fig. [5Jc) without affecting the 
d-mode nature of ip v before and after the bounce. 

Based on the above arguments, we suggest that 'our 
C- and d-mode initial conditions' in Table [TT| produce 
the proper C- and d-modes which are valid throughout 
the bounce, as presented in Figs. QJa), 01b), I3 a ): an d 
[5fb)- In order to obtain these numerical results produc- 
ing proper C- and d-mode decomposition, we have to find 
the precise initial condition; the method will be explained 
in Sec. IIVA3I 

In an ideal situation we may use the initial conditions 
of ip v and <p x in Eq. (|58[) which are valid when the large- 
scale and the adiabatic conditions are satisfied exactly; 
these initial conditions are summarized in Table [TTT1 and 
numerical results will be presented in Sec. IIV A 21 Since 
these two conditions are not exactly satisfied in reality, 
we have to find more accurate initial conditions which 
allow the correct C- and d-mode behavior at the initial 
phase, and preferably throughout the bounce. The cor- 
rect initial conditions we found are presented in Table HT1 
and the numerical results are presented in Figs. H] and [B] 

As noted in Sec. IIII B[ we solve numerically two inde- 
pendent sets of equations based on two different temporal 
gauge conditions: these are Eqs. (|5"Tj) - (j5~4")) in the uniform- 
tr gauge, and Eqs. (j5T)[) - (foT]) in the zero-shear gauge. We 
derive initial conditions for variables in the zero-shear 
gauge using relations between the uniform-cr gauge and 
the zero-shear gauge in Eq. (|62p . Any gauge- invariant 
combination of variables should show the same behav- 
ior independently of the gauge conditions we took. This 
provides a numerical check of the accuracy of our inte- 
gration. Although evolutions of (f v and ip x in Figs. [Hand 
[6] sometimes show singular behaviors near the bounce, 
these are due to the singular behavior of background vari- 



TABLE II: The initial conditions which give the precise (thus, 
correct) C- and d-mode behaviors in both collapsing and ex- 
panding phases away from the bounce. We call these 'our C- 
and dmode initial conditions, or simply 'our C- and d-modes. 
We take k = 0.0005 and the initial 77 = —500. The numerical 
results based on these initial conditions are presented in Figs. 



variable 


C-mode 


d-mode 




-7.8058229814554E- 


-06 


-1.5568456514629E+03 


5<t>5a 


4.6325084877449E- 


■08 


8.9724639867256E+00 


<fi5a 


5.4370235008682E- 


■06 


-8.9884525573535E+02 


Via 


2.7003984444483E- 


■08 


5.1802544967425E+00 


fv 


9.9439752296909E- 


■06 


1.0015746283898E-05 




2.5769163366172E- 


10 


5.8225396948330E-08 



ables multiplied to form these gauge-invariant combina- 
tions, e.g., see Eq. (|4"3T> . In Figs. [7] and [5] we show that 
during our numerical integration the original variables 
used in the integration behave smoothly. 

In Figure [9j we show the C-mode evolution of ip v and 
ip x for a different scale which enter the horizon at later 
epoch; here we set k = 10 -6 whereas in other figures we 
used k = 0.0005, and the initial conditions differ from 
the other figures. 

Our results in Figs. [4] and [6] show that by taking 
the correct initial conditions representing the C- and d- 
modes for ip v and ip x , the identities of C- and d-modes re- 
main valid throughout the bounce. That is, the C-mode 
remains as the C-mode throughout the bounce, and sim- 
ilarly for the d-mode. Although in the collapsing phase it 
is the d-mode which is relatively growing, it decays away 
in the expanding phase. Thus, if we are interested in 
the relatively growing mode in the final expanding phase 
we have to consider the C-mode even in the collapsing 
phase. This is one of the major points we have empha- 
sized in our previous works based on analytic arguments 
in @, H, 0] , and here we confirm it by using the numerical 
integration based on a specific bounce model. 



2. Analytic C- and d-modes 

In Figs. [TOTa) . fTUTb) . [TTT a) . andfTlTb') we present some 
generic behaviors of imprecisely imposed C- and d-mode 
initial conditions. Unless we set up initial conditions pre- 
cisely as in 'our C- and d-modes' presented in Table ITT| 
it is very likely that the curvature perturbations behave 
similar to the ones in Figs . [TOTa) . [TOTb') . [TTT a) . andfTlTb'). 
What is actually presented in these figures is based on 
'analytic C- and d-mode' initial conditions in Table \TU\ 
These initial conditions are based on the large-scale limit 
solutions of adiabatic perturbation in Eqs. (|39[) and (|40[) . 
Thus, these are valid for an ideal situation assuming ex- 
act validity of the adiabatic condition and the large-scale 
limit. 
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FIG. 4: Evolutions of our C-modes for <p v and <p x , and the adiabatic conditions; the large-scale conditions are presented in 
Fig. [5] Panels (a) and (b) show the behaviors of our C-modes for ip v and tp x . We use the initial conditions in Table. [II] Panels 
(c) and (d) show the adiabatic conditions for tp v and tp x , respectively; for tp„ and ip x we plot the LHSs of Eq. (|37[) . Panel 
(c) shows that the adiabatic condition for ip v is broken by sharp peaks (divergences) six times near the bounce; we can show 
that the four sharp peaks at the outer edge and the inner edge surrounding the bounce are caused by z" jz = 0, and the other 
two sharp peaks in the middle are caused by /i + p = 0. The adiabatic condition for tp x in Panel (d) also shows two sharp 
peaks just before the bounce because ip x crosses zero at these two points; see Eq. (137f) and Fig. 04) which shows evolution of 
tp x in a linear scale. Later, in Fig. [T2]we will show that the isocurvature perturbation is not affected by the presence of these 
sharp peaks in the adiabatic condition and in the large-scale condition; the latter will be presented in Fig. [5] Based on this 
observation and the physical behaviors of the curvature variables throughout the bounce, we suggest that such peaks do not 
have physical impact on the C- and d-mode natures of the perturbation. Panel (a) shows a spike in the evolution of cp„ just 
after the bounce; it actually have two spikes at fi + p = 0, see Eq. (|43p and Fig. [7] Away from the bounce we have <p x — (3/5)</?„ 
which is the well known ratio for w = background. 



In Figs. ITUlc), ITOfo), ITTlc). andHUd) w e present the 
adiabatic conditions. The large-scale conditions are the 
same as in Fig. [5] Detailed descriptions of the evolutions 
and the adiabatic conditions can be found in the figure 



captions. Sharp spikes away from the bounce in the adia- 
batic conditions in Figs. fTTJTc) . flUT d) . andQTJd) occur at 
ip v = or ip x — which can be identified in Figs. fTOT a). 
[TUf b) , and [TUb) . The curvature perturbations vanish at 
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(a) sign(ct) ■log(o-/o, min ) (b) sign(ct) -log(a,/a mln ) 

FIG. 5: Large-scale conditions for <p v and ip x . Panels (a) and (b) show the large-scale conditions for ip v and <p x equations; for 
tp v and ip x we plot k 2 /(«" /z), and k 2 /[(1/z)" /(1/z)], respectively. As the large-scale conditions involve only the background 
evolution and the wavenumber, the same large-scale conditions apply to both C- and d-modes independently of the initial 
condition. Panel (a) shows that the large-scale condition for ip v is broken four times by sharp peaks (divergences) near the 
bounce. These divergences occur at z" / z = 0. As explained in Fig. [4] as the isocurvature perturbation is not affected by the 
presence of these sharp peaks in the large-scale condition, such peaks are not supposed to have physical impact on interpreting 
the C- and d-mode natures. 



these points as the dominating mode changes between 
the C-mode and the d-mode. Sharp spikes at the bounce 
in Figs. Unjd) and [TTT d) occur because <p x = at the 
bounce. Even in our d-mode case we have ip x = at the 
bounce, see Fig. [5] Figures [TOTb) and [TTTb) are effec- 
tively the same; although we call Fig. fTUT b) an analytic 
C-mode, apparently, the initial condition is dominated 
by the d-mode from the beginning. 

The adiabatic conditions for (p v infTOFc) and [TlT c) are 
severely broken near the bounce. Evolutions of ip v in 
Figs.QJj^a) and [TIT a) also show switching of the d-mode 
before the bounce into the C-mode after the bounce. Fur- 
thermore, Figures [121 and [T3l show that, for the analytic 
C- and d-modes, the isocurvature perturbation is signif- 
icantly excited at the bounce. Based on these observa- 
tions we conclude that we cannot trace the C-modes af- 
ter the bounce to the d-modes before the bounce in Figs. 
HOT a) and [Til a). This conclusion should be compared 
with the one made in our d-mode case. Although, the 
adiabatic condition for our d-mode is also severely bro- 
ken near the bounce [see Fig.[6]Jc)], there we have argued 
how this could happen for the d-mode of tp v without af- 
fecting the d-mode nature before and after the bounce. 
In our d-mode case, Fig. [SJa) shows that the ci-mode na- 
ture is preserved before and after the bounce, and Fig. [13] 
shows that the isocurvature perturbation is not excited 
near the bounce. 

The evolutions show that each mode is soon dominated 
by the relatively growing mode (the d-mode in a col- 



TABLE III: The analytic initial conditions in Eq. (J58J). This 
initial condition is valid for an ideal situation assuming ex- 
act validity of the adiabatic condition, and taking the lead- 
ing order solution in large-scale expansion presented in Eqs. 
(|39[) and (|40[) . We call these 'analytic C- and d-mode initial 
conditions', or simply 'analytic C- and d-modes'. We take 
k = 0.0005 and the initial n = —500. The numerical results 
based on these initial conditions are presented in Figs. HOI and 

nu 
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lapsing phase, and the C-mode in an expanding phase 
are the relatively growing). In fact, these analytic initial 
conditions correspond to mixtures of our precise C- and 
d-mode initial conditions. Based on these observations 
we conclude that we cannot trace the final C-modes af- 
ter the bounce in all these figures to the ones provided 
initially. 
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FIG. 6: Evolutions of our d-modes for tp v and <p x , and the adiabatic conditions; the large-scale conditions are presented in 
Fig. [5] Panels (a) and (b) show the behaviors of d-modes for ip v and tp x , respectively. We use the initial conditions in Table 
Ull Both ip v and ip x cross zero at the bounce, see Figs. |8jd) and BTe) which show the evolution in linear scale. <p v diverges 
twice near the bounce at n + p = 0, see Eq. (|43[) an d Fig. |8jd). Panels (c) and (d) show the adiabatic conditions for ip v and 
ip x . Although, the adiabatic condition for ip v equation is severely broken near the bounce, we interpret this happens because 
the d-mode of ip v is (k/aH) 2 -order higher in the large-scale expansion compared with the same mode of ip x , i.e, the leading 
order d-mode of ip v has been canceled out; see Eqs. (I39|l and (|40[) . or Eqs. (|41|l and (|42[) . Panel (c) also shows four sharp 
peaks (divergences) near the bounce which are caused by z" / z = 0; compared with Fig. |3Jc), the two sharp peaks potentially 
caused by fi + p = disappeared because the d-mode of (p v diverges in proportion to + e.g., see Eq. (|39[1 and Panel (a). 
The adiabatic condition for tp x has one sharp peak at the bounce because tp x crosses zero at this points; see Eq. (I37|l and Fig. 
|8je). As explained in Fig. |4j the isocurvature perturbation is not affected by the presence of these sharp peaks in the adiabatic 
condition, and the d-mode before the bounce remains as the d-mode after the bounce. Thus, such peaks do not have physical 
impact on the d-mode natures of the perturbations. 



3. Method of finding our C- and d-modes 

Here, we explain how we found 'our C- and d-mode 
initial conditions' which lead to correct behaviors of the 



C- and d-mode throughout the bounce. Since the d-mode 
is rapidly growing in a collapsing phase, thus dominat- 
ing, we can easily find d-mode in that phase; similarly 
it is easy to find C-mode in an expanding phase. How- 
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ever, unless we impose a quite precise initial condition, it 
is very likely that, the relatively growing mode in that 
phase soon dominates the evolution. Thus, after the 
bounce a d-mode before the bounce is likely to switch 
to a dominating mode in an expanding phase which is 
the C-mode. We have shown that our C- and d-modes 
can be identified as the proper C- and d-mode which pre- 
serve their nature before and after the bounce. An im- 
precise initial condition, which tends to reproduce only 
the dominating solutions, can be regarded as a mixture 
of our (precise) C- and d-mode initial conditions. In Sec. 
IIV A 21 we showed that such a mixed initial condition 
leads to the break down of the adiabatic condition near 
the bounce and excitation of the iscurvature perturba- 
tion. Without a precise initial condition, it is hard to 
maintain d-mode nature after the bounce. Similarly, it 



is quite difficult to find C-mode initial condition which 
maintains its C-mode nature during the collapsing phase. 

In order to find the precise initial condition for the 
d-mode we begin with the analytic initial condition im- 
posed close to the bounce, and integrate forward and 
backward in time; the analytic initial condition is pre- 
sented in Eq. ([551) which is based on Eqs. and (gDJ). 
Our d-mode initial condition found in this way is pre- 
sented in Table [TTJ As we have shown this d-mode initial 
condition leads to "the d-mode" which maintains its d- 
mode natures even after the bounce; see Fig. El^a) and 
|6jb). Since the analytic initial condition in Eq. (l58l) 
assumes the exact adiabatic condition, the isocurvature 
perturbation of our d-mode naturally vanishes near the 
bounce, i.e., — 0; see Fig. ITcH 

Our C-mode initial condition is also found similarly. 
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FIG. 7: The raw behaviors of gauge-invariant variables numerically solved using Eqs. (|51[) -(|54 [ ) f° r C-modes of tp v and ip x . 
Panels (a)-(c),(e) show that the variables we used in numerical integration behave smoothly throughout the evolution. We 
obtain the behaviors of tp v and <p x in Panels (d) and (f) using Eq. (|43p . Panel (c) magnifies Panel (b) near the bounce region. 
This figure explains why we have only one sharp spike after the bounce visible in Panel (d) and Fig. HJa); the spikes are caused 
by vanishing fj, + p in Eq. (g3J. The arrows indicate the epochs where fi + p — 0. 



The C-mode behavior of <p v can be found if we begin with 
the analytic initial condition imposed close to the bounce, 
and integrate forward and backward in time. The be- 
havior of ip x found in this way, however, shows d-mode 
behavior shortly before and shortly after the bounce. In 
order to find the C-mode behavior of ip x without ap- 
pearance of the d-mode behavior near the bounce as in 
Fig. SJb) , we should begin with the analytic initial con- 
dition imposed away from bounce, and integrate forward 
and backward in time. In order to have the C-mode be- 
havior in Fig. 01 we used the initial condition found by 
imposing the analytic initial condition at a point where 
S(j, a vanishes (we call it an = condition) in Fig. 
[T^c). In this way, behaviors of C-modes in Fig. [5] how- 
ever, depend on whether we found our initial condition 
by imposing S^, a — before or after the bounce as in 
Figs. frW c) and fT4T f). The C-mode behaviors in Figs. [4] 
and [7] are related to Fig. H^c). In Fig. [14] we present 
the counterparts of Figs. H] and [7] now related to Fig. 
[T4jf). Comparison of these figures will show that near 
the bounce Fig. Q3] is the reflection of Figs. 0] and [7] at 
the bounce. In Fig. rjS] we present an additional case 
where we begin by imposing the analytic initial condi- 
tion at one of the epochs of vanishing /i + p. In this way, 
tf x — at only one point, see Fig.rjS^b) andfTHe). If we 
begin by imposing Sfa = condition closer to or more 
away from the bounce we will have d-mode behavior of 
ip x near the bounce region. Thus, the C-mode behaviors 
in Figs. 0Jb), HHb), and fTSTb) are closest to the analytic 
behavior of C-mode of ip x in Eq. ([30]) ■ 



B. Tensor- type perturbations 

The evolution of tensor-type perturbation is described 
by Eq. ([64]) . Figures fTBlTLSl show diverse behaviors of 
the tensor-type perturbation depending on initial condi- 
tions. Using the C Q( a-mode initial condition in Eq. (|65j) . 
Fig. [16] shows a conserved evolution of C^ throughout 
the bounce: thus, we can clearly identify this as the C' a p- 
mode in Eq. (|65[) . However, although we naively expected 
to be able to find a precise d a( a-mode initial condition in 
the collapsing phase which will be preserved as the same 
d a p-m.o&e even in the expanding phase, the numerical re- 
sult in Fig. 1 171 shows that the d a ^-mode initial condition 
switches to the C a( 3-mode. 

This result differs from our following naive anticipa- 
tion. The equation for tensor-type perturbation is quite 
similar to the one for tp v in an adiabatic limit: compare 
Eq. ([M]) with Eq. ([50]) . In large-scale limits we have an- 
alytic solutions in Eqs. ([65]) and ([39]) . In case of d-mode 
for the scalar-type perturbation in Fig. [SJa) we presented 
a proper c?-mode evolution before and after the bounce. 
Thus, previously we anticipated that even for the tensor- 
type perturbation similarly we can find the proper ini- 
tial condition which produces the d a( g-mode before the 
bounce continuing as the d a p-mode after the bounce. In 
the following paragraph we can show that it is not possi- 
ble to have such a behavior. Thus, any d a p-mode initial 
condition in a collapsing phase is later dominated by an 
apparent C Q( a-mode solution in an expanding era. This 
implies that we do not have the proper d^-mode which 
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preserves its nature throughout the bounce. Figurc fTTl b) 
shows that the initially supplied d^-mode still survives 
as the same <i Q ^-mode even in the expanding era. Figure 

[P7T a) shows that, after the bounce, C^i is simply dom- 
inated by an additional Cc^-mode appearing after the 
bounce; this apparently comes from the lower-bound of 
integration of Eq. ([65]) . 

Here we can show why it is not possible to find the 
initial condition which remains as the d Q( g-mode after 
the bounce. // the e? Q ^-mode in the collapsing phase 
could somehow switch to the same mode in the expanding 
phase, at the bounce C^l' should either change its sign 
or cease to be differentiable. However, from Eq. |65|) . for 
k = 0, we have 



C 



(ty 

af3 



d a/3 (k) 



(70) 



where d a p is a temporal constant. Thus, once the sign 
of d a {3 is decided, then its sign is fixed throughout the 

evolution, and cj£l' is always differentiable. In case of 
the large-scale solution of ip v in Eq. ([39|) . changing sign 
of z 2 plays an important role in explaining the different 
behaviors between d-mode of ip v and J Q /3-mode of C^; 
Fig. [5Jd) shows that <p v does not only change its sign 
but also is not differentiable in a couple of places. For a 
non-vanishing k, from Eq. (|66p we have 
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(71) 



Since we are considering the d a p-mo&e solution in the 
large-scale limit, the large-scale second-order correction 
term (the second term in RHS) in Eq. (|7ip should be 
negligible compared with the leading order term in Eq. 



20 




(d) 



-2 2 

I 

sign(a>) ■log(a,/a, min ) 



4000 



2000 



-2000 



-4000 





I i ii 
r ~ 

I 1 


, , , 


i i I,, 



(e) 



sign(a>) •log(ot//a> 



FIG. 8: The same as Fig. [7] now, for the d-mode. Panels (a)-(c) show that the variables we used in numerical integration 
behave smoothly throughout the evolution. 



TABLE IV: The initial conditions for tensor-type perturba- 
tion with k = 0. For k = the large-scale condition is met 
exactly. We take initial rj = —500. Results based on these 
initial conditions are presented in Figs. 1161171 
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([70]) : thus, the sign of C^l' cannot be changed in that 
limit. This confirms that our numerical result for the 
d a /3-m.ode presented in Fig. I17f a) is unavoidable. 

In Fig. 1181 we present an evolution of initially near 
Ca/3-mode which transits to the d a /3-mode after the 
bounce. Compared with the precise (7^-mode in Fig. 
1161 in this Figure we introduced a small amount of ini- 
tially negligible but non-vanishing d a p-mode. The time 
derivative of C£p in Fig. IMb) shows that the apparent 

do^-mode in the expanding phase is exactly the one pro- 
vided initially. After the bounce, the initially dominating 
Ca/3-mode has been canceled out by a contribution from 
the lower-bound of integration of Eq. (f6"5|) . 



V. POWER SPECTRA 

A. Scalar-type perturbations 

As the origin of seed fluctuations for the large-scale 
structure we consider quantum fluctuations of the fields 



and the space-time metric ever present in the collapsing 
phase. Assuming certain dynamics of the scale factor 
we have the perturbation equations described by Bessel 
equations with the solutions in Hankcl functions as in 
Eq. (|49p . The canonical quantization process well known 
in the literature leads to the same equations and solu- 
tions applicable to the quantum fluctuations. In order 
to correspond the large-scale solutions with the C- and 
d-mode decomposition in Eqs. (f3T))) and PO)) to the exact 
solutions with Hankel functions in Eq. (|49]) , we follow 0] . 
The power spectrum and spectral index are defined as 
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(72) 



where x can be ip v (k 1 r]) or <p x (k, rj) in the Fourier space. 

In the power-law case away from the bounce where the 
ghost field is negligible, Eqs. ([50]) and pTj) lead to Bessel 
equations. Using the quantization based on the action 
formulation, we have the exact mode function solutions 
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FIG. 9: Long term evolution of the C-mode. This figure shows 
the evolutions of C-modes using k = 10 -6 . The inset figure 
shows that the scale becomes subhorizon near log (a/a m i n ) 
1 0. Both Lp v and tp x begin oscillations as the scale comes 
inside the (visual) horizon. In a realistic situation, the horizon 
crossing of relevant cosmic scales should occur long after the 
bounce. 



The Hankel functions can be expanded as [2z 



H^\k\n\) =J2-< 



(k\n\? 





(74) 



In large scale limit, the first (second) term in the paren- 
thesis dominates for v < (y > 0). In Eqs. (|39|) and 
([40]) the leading orders of the C-modes are time indepen- 
dent whereas the leading orders of the d-modes behave 
as <p v oc l^l 2 " 1 ' and ip x oc \?]\ 2l/u . Since 



(75) 



we can easily identify the first and the second terms in the 
parenthesis of Eq. ([74]) as the d-mode and the C-mode, 
respectively. 

After taking the simple vacuum choice C\ = \ and C2 = 
0, the spectral indices for the C-mode and d-mode of tp v 



and ip x can be read as 
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(76) 



Notice that for the C-mode, the spectral index of tp v and 
ip x coincide. 

The near Harrison-Zel'dovich scale-invariant spectrum 
corresponds to ns — 1 — which is consistent with cur- 
rent observations of the cosmic microwave background 
radiation [23| and the large-scale structures [24(. The 
ekpyrotic or the cyclic scenario has 0<(j<l, thus from 
Eq. (|76p these scenarios have a scale-invariant spectrum 
for d-mode of tp x , but these have a quite blue spectrum 
with ns — 1 — 2 for the C-mode [8] ; incidentally, we have 
ns — 1 — 4 for the d-mode of tp v . Although the d-mode 
of ip x shows a scale-invariant spectrum, we have shown 
that the G?-mode eventually decays away in the expanding 
phase. Thus, we cannot see this spectrum in the current 
expanding phase of the universe. The observationally rel- 
evant one is the C-mode spectrum. Thus, for the ekpy- 
rotic or the cyclic scenarios we have ns — 1 — 2 which is 
too blue compared with observations. Here, we assume 
a smooth nonsingular bounce with valid linear pertur- 
bation through the bounce; this differs from the original 
ekpyrotic/cyclic scenarios where the effective four dimen- 
sional spacetime goes through a singularity In the 
pre-Big Bang scenario we have q = 1/3, thus ns — 1 = 3 
for the C-mode Q. 

In our simple bounce model case with q = 2/3, we 
have a scale-invariant spectrum for d-mode of ip v , but 
the d-mode is meaningless because it decays away in an 
expanding phase. Our model shows a quite blue spectra 
for C-mode with ns ~ 1 = 6. In Figs. [1^112^1 we present 
the power spectra of C- and d-modes of ip v and <p x based 
on numerical integrations; we also present characteristic 
evolutions of ip v and lo x . The initial spectra are based on 
the large-scale limit of Hankel function solutions in Eq. 
(|73p which have their origins in quantum fluctuations in 
the collapsing phase. The spectral indices of final power 
spectra coincide with our analytic estimations made in 
Eq. ([7S]>. In Figs. EMM we present the initial and fi- 
nal power spectra from the analytic C- and d-modes to- 
gether, and the characteristic evolutions of (p v and tp x . 

Based on numerical evolutions of certain initial condi- 
tions which are similar to the 'analytic C- and d-modes' 
of ours, the authors of 12| reported that the d-mode in 
the collapsing phase survives as the C-mode in the ex- 
panding phase, see Fig. Illf a); similar results were also 
reported in @, [l3|, [3] • This result leads these authors to 
claim that for q = 2/3 the quantum fluctuations in the 
collapsing phase give rise to a scale- invariant spectrum in 
an expanding phase. Our numerical result in Fig. l2"TT c) 
apparently confirms this result. However, based on the 
behaviors of our C- and d-modes, the initial conditions 
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FIG. 10: Evolutions under the 'analytic C-mode' initial condition in Table Hill which is valid for an ideal situation based on an 
analytic solution. This figures is the same as Fig. |4j now for an analytic C-mode initial condition. Panels (a) and (b) show 
evolutions of tp v and ip x . Panels (c) and (d) show adiabatic conditions for ip v and tp x . The large-scale conditions are the same 
as in Fig. [5] Compared with 'our C-mode' presented in Fig. [4] which shows that the C-mode nature is maintained before and 
after the bounce, the present 'analytic C-mode' leads to mixed behaviors of the C- and d-modes. In the present case, the initial 
condition in Table Hill corresponds to a mixture of our precise C- and d-modes presented in Figs. [4] and [6] respectively. Panels 
(c) and (d) show that the present initial conditions lead to failures of the adiabatic conditions both for tp v and <p x during the 
evolution. Two sharp spikes away from the bounce in Panels (c) and (d) are caused because the curvature variables <p v and 
tp x vanish as the dominating mode switches between the C-mode and the d-mode, see Panels (a) and (b). The sharp spike at 
the bounce in Panel (d) is caused because ip x vanishes at the bounce; apparently, the behavior of ip x shows its d-mode nature 
before and after the bounce, thus vanishes at H = 0, see Eq. (|40|l . As in our C-mode in Fig. [4] such a sharp peak does not 
affect the d-mode nature of ip x at the bounce. The adiabatic condition for tp v in Panel (c), however, is severely broken near 
the bounce. The evolution of <p v in Panel (a) also shows switching of the d-mode before the bounce into the C-mode after the 
bounce. Furthermore, later in Fig.[l2]we will show that, for the analytic C-mode, the isocurvature perturbation is significantly 
excited at the bounce. Based on these observations we conclude that we cannot trace the C-mode after the bounce to the one 
provided initially. 
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FIG. 11: The same figures as in Fig. 1101 now for the analytic d-mode initial condition in Table Hill Compared with 'our d- 
mode' presented in Fig. [6] which show that the d-mode nature is maintained before and after the bounce, the present 'analytic 
d-mode' leads to mixed behaviors of the C- and d- modes. In the present case, the initial condition in Table Mil corresponds 
to a mixture of our precise C- and d-modes presented in Figs. [4] and [6] respectively. Panels (c) and (d) show that the present 
initial conditions lead to failures of the adiabatic conditions both for ip v and <p x during the evolution. One sharp spike away 
from the bounce in Panel (d) is caused because ip x vanishes as the dominating mode switches from d-mode to C-mode, see 
Panel (b). As explained in Fig. 1101 the sharp spike at the bounce in Panel (d), caused by vanishing <p x at H = 0, does not 
affect the d-mode nature of ip x at the bounce. As in Fig. 1101 the adiabatic condition for tp v in Panel (c) is severely broken near 
the bounce. The evolution of ip v in Panel (a) also shows switching of the d-mode before the bounce into the C-mode after the 
bounce. Later in Fig. [T3]we will show that, for the analytic d-mode, the isocurvature perturbation is significantly excited at 
the bounce. Based on these observations we conclude that we cannot trace the d-mode after the bounce to the one provided 
initially. The adiabatic condition for ip v is also severely broken for 'our d-mode' in Fig. [S] In the case of our d-mode, in Fig. |S] 
we have argued that this happens because the d-mode of ip v is (k/aH) 2 -order smaller than the d-mode of ip x , thus quite small, 
see Eq. (|39[1 . Furthermore, in the case of our d-mode the isocurvature perturbation is not excited near the bounce, despite the 
apparent breakdown of the adiabatic condition for cp v , see Fig. 1131 
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used in these analytic C- and d-modes are mixture of the 
pure (our) C- and d-modes. Under these initial condi- 
tions, as the evolution proceeds we have shown that the 
isocurvature perturbation is also excited near the bounce, 
see Fig. [13] Thus, it is not correct to interpret that the 
final C-mode spectrum in the expanding phase is gener- 
ated from the initial d-mode spectrum in the collapsing 
phase. In our interpretation, although the evolution be- 
gins with near d-mode (which is actually a mixture of 
our, thus precise, C- and d-modes), as we approach the 
bounce not only the C- and d-modes but also additional 
two isocurvature-type perturbations are all excited; this 
eventually leaves only the dominating mode (which is 
the C-mode in an expanding phase) after the bounce. 
In any case, our numerical result in Fig. I21f c) shows 
that near Harrison-Zel'dovich spectrum can be generated 
from the analytic d-mode initial condition with the help 
of the isocurvature perturbation simultaneously excited 
near the bounce. In case of our C- and d-modes where 
the C- and d-mode natures are preserved before and af- 
ter the bounce, and the isocurvature perturbation is not 
excited. In that case, near Harrison-Zel'dovich spectrum 
in the d-mode of <p v is observationally not relevant. 

A scale-invariant spectrum for the C-mode can be ob- 
tained for q 1 which is the power-law inflation in an 
expanding phase. In the collapsing phase, however, q > 1 
implies that the relevant scales come inside the horizon 
near the bounce, thus invalidating the large-scale condi- 
tion we used to get the above result. In the sub- horizon 
scale the modes oscillate, and lose their previous large- 
scale identity of the C- and the d-modes [8j. From Eq. 
(|76p . for any model with q < 1 during the quantum gen- 
eration stage, we have n$ — 1 > 0, thus having a blue 
spectrum. 

Therefore, we conclude that while the large-scale con- 
dition is satisfied and the adiabatic condition is met 
(i.e., the isocurvature perturbation is not excited) during 



the bounce, it is not possible to obtain near Harrison- 
Zel'dovich scale-invariant density spectrum through a 
bouncing world model as long as the seed fluctuations 
were generated from quantum fluctuations of the cur- 
vature perturbation in the collapsing phase. This con- 
clusion assumes that Einstein gravity is dominating 
throughout the bounce including the quantum genera- 
tion stage. 



B. Tensor- type perturbations 

The quantum fluctuations of the space-time metric 
detached from the fields provide seed fluctuations for 
the tensor- type perturbation (gravitational waves). For 
tensor-type perturbation, the power spectrum and spec- 
tral index are defined as 
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where we introduced a decomposition based on the two 
polarization states (2f| 

C$(x,t) ^ WolJ ^0c3(x,t;k) 

S ^/(^? eikX " ik(<)eS(k) ' 
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FIG. 12: Evolutions ol the isocurvature perturbation variable for 'our C-mode' (solid line) and the 'analytic C-mode' (dotted 
line). Panel (b) magnifies Panel (a) near the bounce. Panel (c) shows Panel (a) in logarithmic scale. The sharp divergences at 
the bounce occur because we have <f> = at the bounce. For our C-mode, S$ a nearly vanishes near log (a/a m m) — — -2; this 
is because we have found our C-mode initial condition beginning at this epoch, see the text for explanation. Although Figs. 
Uc) and (d) show that the adiabatic conditions for our C-modes of ip v and <p x are broken by sharp peaks near the bounce, 
the above figures show that the isocurvature variable show smooth behavior near the bounce. As explained in Fig. [4] such 
sharp divergences occur because we have fi + p = near the bounce and ip x crosses zero twice near the bounce, see Fig.^b). In 
case of the analytic C-mode, Figs. UOf c) and (d) show that the adiabatic conditions are severely broken near the bounce, and 
are sharply broken in two places away from the bounce. As explained in Fig. 1101 the two sharp spikes away from the bounce 
occur because of the change in dominating mode between C-mode and d-mode, where the curvature variables ip v or ip x vanish. 
The above figures show that for the analytic C-mode the isocurvature variable significantly changes its behavior at the bounce. 
Based on these behaviors, we argue that for our C-mode initial condition although the adiabatic conditions are broken by 
sharp peaks near the bounce, because the isocurvature perturbation shows the smooth behavior, the isocurvature perturbation 
does not give physical impact on the adiabatic nature of our C-mode throughout the bounce. Thus, although our C-mode 
initial condition shows apparently non- vanishing isocurvature mode (see Table 1X1)1 . we argue that our C-mode initial condition 
keeps the adiabatic nature throughout the bounce. On the other hand, for the analytic C-mode, the adiabatic conditions are 
severely broken near the bounce, and the isocurvature perturbation shows sharp change at the bounce. Thus, we argue that 
the isocurvature perturbation gives physical impact on the adiabatic nature of the analytic C-mode; this analytic C-mode is 
based on an initial condition with 'precise' adiabatic condition (see Table ITTT)) . That is, for the analytic C-mode, although 
the initial condition is imposed assuming the adiabatic condition is exactly met, as the evolution proceeds the isocurvature 
perturbation is also excited, especially near the bounce. Thus, the numerical evolution shows that after the bounce the analytic 
initial C-mode in fact becomes a substantial mixture of the adiabatic and isocurvature modes. 
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In order to correspond our large-scale solution in Eq. 
(1551) to the exact solution in Eq. Q68p. we follow a similar 
reasoning used for the <p v . In Eq. (frj5"| the leading order 
of the C a p-mo&e is time independent whereas the leading 
order of J^-mode behaves as cfp oc \rj\~ 2Vg . Since 



OM) °c hr«<' 2 ^M), (79) 

we can easily identify the first term in the parenthesis of 
Eq. ([74]) as the C^-mode and second term as the d a p- 
mode. We can read the spectral indices for C a p and d a p 



as 



n T,C a 
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(80) 



These spectral indices of tensor-type perturbation coin- 
cide with the spectral indices of the ip v in Eq. (|T6[) . 

From the classical evolution of the tensor-type pertur- 
bation studied in Sec. lIVBl we notice that both the initial 
Cap- and d a/ 3-modes survive as the C^-mode after the 
bounce. If both modes are generated from the quantum 
fluctuations with comparable amplitude, as the collaps- 
ing phase proceeds the initial C a p-mode becomes negligi- 
ble compared with the initial <i a ^-mode. Thus, after the 
bounce the dominant C^-mode should have origin from 
the initial d ai g-mode; as the d^-mode grows rapidly in 
the collapsing phase, its final amplitude surely depends 



2G 




on the duration of the collapse after its quantum origin. 
This is in contrast to the situation of the scalar-type per- 
turbation where the initial d-mode remains as the d-mode 
even after the bounce, thus decays away in the expanding 
phase. 

In the pre-big bang scenario with q = 1/3, we have 
n T,c al 3 = ^ and n T j a/3 = 3. In the ekpyrotic or cyclic 
models with q ~ 0, we have ji r g o(J = 2 and n T = 4. 

Notice that it is the d a( 3-mode initial condition in the 
collapsing phase which survives as the dominating C a p- 
mode after the bounce in case of the tensor-type pertur- 
bation. Thus, the dominating C^g-mode perturbation 
in expanding phase produced in the ekpyrotic and cyclic 
scenarios has tit = 4 spectrum; this also differs from 
riT = 2 which is claimed by the authors of these sce- 
narios [H|. Meanwhile, in our simple bounce model with 
q = 2/3, we have = 6 and n Ti g = 0. From Eq. 

(|80| we notice that only for q — 2/3, thus w = 0, we have 
a scale-invariant near Harrison-ZePdovich type spectrum 
n T g = for the tensor-type perturbation. In Fig. [23] 
we present the final spectra of the initial C a p- and dap- 
modes in case of w = 0; these numerical results coincide 
with our analytic estimation made above. 



VI. DISCUSSION 

In order to explain the origin of large-scale structure in 
the universe, several models of bouncing scenarios 0, 0| 
were suggested as possible alternatives to the inflation. 
The pre-big bang scenario in 0] failed to produce the ob- 
servationally required scale-invariant spectrum for scalar- 
type perturbation; this model has q = 1/3 thus leads to 
(ris — l)tp v ,c — 3 and n T ^ ~ 3. Whereas, the authors 
of ekpyrotic or cyclic scenarios [H claimed that their sce- 
narios successfully produce the required spectrum. The 



latter authors find that the quantum fluctuations of d- 
mode of tp x show near Harrison-Zel'dovich scale- invariant 
spectrum and claimed that this mode gives rise to the 
C-mode after the bounce. As long as the adiabatic con- 
dition and the large-scale condition are met during the 
bounce (and, of course, assuming the linear perturbation 
theory survives the bounce) from the analytic point of 
view we can see that d-mode before the bounce remains 
as the same d-mode after the bounce and similarly for the 
C-mode. Since the d-mode is transient in an expanding 
phase, we have to consider the C-mode as the proper seed 
for the large-scale structure. For the ekpyrotic or cyclic 
scenarios with < g « 1 the quantum fluctuations of 
the C-mode of both ip x and tp v show very blue spectra 
with ns — 1 — 2, thus fail to explain the observation. 

In this work we have confirmed the above expecta- 
tion based on analytic arguments made in Q by follow- 
ing the C- and d-modes numerically in a simple non- 
singular bouncing world model. We found the C- and 
d-modc initial conditions in which the nature of the 
C- and d-modes are preserved through our nonsingular 
bounce world model. Thus, we have shown that the scale- 
invariant spectrum claimed by the proponents of ekpy- 
rotic or cyclic models is not possible if the bounce occurs 
non-singularly as in our toy bounce model. In the orig- 
inal ekpyrotic/cyclic scenarios based on five dimension 
the effective four dimensional spacetime goes through a 
singularity [H , and the situation cannot be handled using 
perturbation theory in Einstein gravity. Our results also 
imply that it is not possible to generate near Harrison- 
Zel'dovich scale-invariant spectrum from the quantum 
fluctuations in the collapsing phase as long as the isocur- 
vature perturbation is not excited significantly near the 
bounce. We state again the assumptions we made to 
get such a conclusion: (i) Einstein's gravity and linear 
perturbation theory work throughout the evolution, (ii) 
quantum fluctuations of a minimally coupled scalar field 
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FIG. 13: Evolutions ol the isocurvature perturbation variable S^ a for our d-mode (solid line) and the analytic d-mode (dotted 
line). Panel (b) magnifies Panel (a) near the bounce. Panel (c) shows Panel (a) in logarithmic scale. The sharp divergences 
at the bounce occur because we have cj> — at the bounce. For our d-mode S^a nearly vanishes just before bounce near 
log (a/a m i n ) ~ 0.; this is because we have found our d-mode initial condition beginning at this epoch, see the text for explanation. 
Figure [S](c) shows that the adiabatic conditions for our d-modes of ip v are severely broken near the bounce. However, in Fig. 
[6]we argued that the d mode of ip v is (k/aH) 2 -order higher in the large-scale expansion compared with the d-mode of ip x ; i.e., 
the leading order d-mode of tp v has been canceled out, see Eqs. (|39p - (|42p . Notice that in Fig. [Sfd) the adiabatic condition for 
d-mode of ip x is well satisfied. Based on this observation, in Fig. [6] we have argued that despite the apparent severe breaking 
of the adiabatic condition for ip v the adiabatic nature of our d-mode initial condition is preserved throughout the bounce. In 
order to reinforce our such an argument in this figure we show the evolution of isocurvature perturbation variable S^u. The 
above figures show smooth behavior of our d-mode near the bounce. In case of the analytic d-mode, Figs. [TT] (c) and (d) show 
that the adiabatic condition for tp v is severely broken near the bounce, and the adiabatic condition for ip x is sharply broken at 
the bounce and away from the bounce. As explained in Fig. 1111 one sharp spike for ip x away from the bounce occurs because of 
the change of dominating mode between d and C mode, where the curvature variable ip x vanishes. However, in contrast with 
our d-mode, the above figure shows that for the analytic d-mode the isocurvature perturbation variable significantly changes 
its behavior at the bounce. Based on these behaviors we argue the following. Although the adiabatic condition for ip v of our 
d-mode initial condition is severely broken near the bounce, because the isocurvature perturbation shows the smooth behavior, 
the isocurvature perturbation does not give physical impact on the adiabatic nature of our d-mode throughout the bounce. 
Thus, although our d-mode initial condition shows apparently non- vanishing isocurvature mode (see Table HI)) . we argue that 
our d-mode initial condition keeps the adiabatic nature throughout the bounce. On the other hand, for the analytic d-mode 
both the adiabatic conditions are severely broken near the bounce and the isocurvature perturbation shows sharp change at the 
bounce. Thus, we argue that the isocurvature perturbation gives the physical impact on the adiabatic nature of the analytic 
d-mode; this analytic d-mode is based on an initial condition with 'precise' adiabatic condition (see Table ITTT)| . That is, for the 
analytic d-mode, although the initial condition is imposed assuming the adiabatic condition is exactly met, as the evolution 
proceeds the isocurvature perturbation is also excited, especially near the bounce. Thus, the numerical evolution shows that 
after the bounce the analytic initial d-mode in fact becomes a substantial mixture of the adiabatic and isocurvature modes. 



in a collapsing phase provide seed curvature (adiabatic) 
mode fluctuations for the large-scale structure, (iii) the 
large-scale conditions are met during the bounce, (iv) 
the adiabatic conditions are met during the bounce. We 
have carefully examined that the latter two conditions are 
often violated subtly near the bounce without affecting 
curvature perturbation natures in the large-scale limit. 

In our simple bounce model we introduced a ghost 
field a only to have a natural bounce in the background. 
The presence of additional ghost field in addition to the 
ordinary minimally coupled scalar field <fi (which is in- 
troduced to have curvature- type seed fluctuations), how- 



ever, forces us to handle a fourth-order differential equa- 
tion for the scalar-type perturbation. In order to have 
a nonsingular bounce we need a to be important only 
near the bounce. However, the presence of ghost field 
could potentially generates the isocurvature-type pertur- 
bation S^ a due to the perturbed ghost field near the 
bounce. Our above conclusion forbidding near Harrison- 
Zcl'dovich spectrum from the bouncing model is based 
on assuming no isocurvature perturbation being excited 
near the bounce (i.e., the adiabatic conditions are well 
met during the bounce). Under this condition we have 
found the C- and d-modes which preserve their nature 
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before and after the bounce. 

By using the analytic C- and d-mode initial conditions 
we have shown that it is actually possible to have near 
Harrison-Zel'dovich spectrum from the analytic d-mode 
initial condition of ip v for q — 2/3, see Fi gs. I2"l7c) and 



2TTd); this result was reported in [6|, |12|, |13j, Il4j |. In 
this case, although Fig. l2~TTd) shows that the initial d- 
mode before the bounce is switched into an apparently 
C-mode after the bounce, such a transition occurs be- 
cause isocurvature perturbation is significantly excited 
near the bounce; i.e., under such an initial condition all 
four modes (two C- and d- curvature modes and two 
isocurvature modes) are excited near the bounce, thus 
mixing modes before and after the bounce. Our numer- 
ical results show that, unless we impose a rather precise 
initial conditions, it is very likely that we can easily have 
only relatively dominating mode during the evolution: 



these are the d-mode in a collapsing phase, and the C- 
mode in an expanding phase. Figures flWa) . IMb) . mTa) . 
and mT b) based on the analytic C- and d-modes clearly 
show such mixed behaviors between the C- and d-modes. 
We have shown that these are because the analytic C- 
and d-mode initial conditions can be viewed as mixtures 
of our (precise) C- and d-modes, i.e., imprecise. In ad- 
dition, even the isocurvature perturbation caused by the 
presence of the ghost field is excited near the bounce. 

We have shown that for the tensor-type perturba- 
tion, however, it is the d^-mode initial condition which 
survives as the dominating C Q( g-mode in the expanding 
phase. Thus, for the ekpyrotic or cyclic scenarios with 
< q <C 1, we have n T l j afj — 4; this also differs from the 
predictions made by the original authors of these scenar- 
ios who claimed the dominating mode. 

In conclusion, as long as the large-scale condition and 
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(e) sign(a,) ■log(a//a/ min ) (f) sign(a,) ■log(a,/a/ min ) 

FIG. 14: The same C-mode as in Figs.|4]and[7]now using slightly different way of finding our C-mode initial condition explained 
in the text. Panels (a)-(d) correspond to Figs. |4ja)-(d). Panel (e) corresponds to Fig.[7Jf). Panel (f) corresponds to Fig. I12f c') . 
Near the bounce these figures show reflection symmetry at the bounce compared with Figs. [3] [TJf ) , and ll2f c). Figs. [TJa) , [7£b) , 
and[7Je) almost remain the same. 



the adiabatic condition are satisfied, we numerically show 
that the C-mode in a collapsing phase leads to the ob- 
servationally relevant quantity in the expanding phase. 
The d-mode, in contrast, remains as the same d-mode 
even after the bounce. Thus, under these conditions, it 
is not possible to obtain a scale-invariant spectrum in the 
bounce cosmological model based on the quantum fluc- 
tuations in the collapsing phase. Equation (|76|) shows 
that only for 5 » 1 during quantum generation stage 
the bounce model can produce near Harrison-Zel'dovich 
scale-invariant density spectrum. However, since this 
leads to a violation of the large-scale condition near the 
bounce we cannot use this result; the modes are mixed 
up and begin to oscillate, thus we cannot trace the C- 
and d-mode nature through the bounce. That is, as the 
scales come inside the horizon, perturbations start oscil- 
lations and lose their large-scale memory of the original 
C- and d-modes. For the tensor-type perturbation, it is 
the initial d Q/ 3-mode which survives as the dominating 
C a( 3-mode in the expanding phase. 

We may summarize the final spectral indices for the 
scalar- and tensor-type perturbations as 

2 4 - 6q . , 

(n s - 1 )c = , n T d = , (81) 

1 — q ' ap 1 — q 

where we have a oc \t\ q during the quantum generation 
stage before bounce. Meanwhile, if the initial seed fluctu- 
ations were generated from quantum fluctuations during 
expanding phase (without or after bounce), the observa- 
tionally relevant perturbations have 

2 

(n s - l) c = — — - =n Tt c a0 , (82) 



which is the well known result in the literature 0, l26| . 
Notice the difference in the gravitational wave power 
spectrum generated between the bounce and the infla- 
tion scenarios. 

Our analysis in this work is based on a specific non- 
singular bounce model with two scalar fields. Using this 
specific model we were able to follow numerically the evo- 
lutions of scalar- and tensor-type perturbations through 
the bounce. We have shown that the nature of large- 
scale decomposition of the adiabatic solutions far away 
from the bounce, i.e., the C- and d-modes, is preserved 
throughout the bounce as long as the two conditions (the 
large-scale and adiabatic ones) are met. Throughout the 
bounce we have numerically monitored that these condi- 
tions are physically met as long as we take the correct 
and precise initial conditions which assure the C- and d- 
mode natures; we also have monitored the evolution of 
isocurvature perturbation. In order to help readers to 
reproduce our numerical results, in Table HT1 we have pro- 
vided the precise initial conditions we used. Because the 
above argument was originally based on the analytic so- 
lutions with general equation of state or field potential, 
our numerical analysis based on specific bounce model 
may have more generic implications which go beyond the 
specific model we used. Our concrete numerical study 
assures that ourprevious conclusions based on analytic 
reasons in 0, H, are valid. 

In this work, we have assumed Einstein's gravity 
and spatially homogeneous-isotropic background world 
model with linear perturbations are valid throughout 
the bounce. In realistic situation, however, it is likely 
that in the collapsing phase as the model approaches the 
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bounce the perturbations of all three types grow rapidly 
and eventually become nonlinear, thus invalidating the 
spatially homogeneous-isotropic Friedmann background 
world model. Furthermore, as the scale factor shrinks 
we reach the high energy scale where quantum correc- 
tions of Einstein's gravity likely lead to serious correc- 
tion terms in the gravity itself. Although we have used a 
ghost field as the bounce agent in Einstein's gravity, it is 
also likely that the quantum correction terms can cause 
nonsingular bounce more naturally [TlT ] . In this work we 
have assumed that the linear perturbation theory sur- 



vives throughout the bounce which certainly restricts the 
allowed bounce model. 
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FIG. 16: Evolution of C a /3-mode tensor-type perturbation. 
The behavior of C a /3-mode is similar to the C-mode of tp v . 
The Cc/3-mode nature of this mode is preserved throughout 
the bounce; i.e., the time derivative of vanishes through- 
out the evolution. 
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FIG. 17: Panel (a) shows the behavior of initial d a p-mode which apparently transits to the C Q ^-mode after the bounce. It is 
not possible to find the initial condition which maintains the d a p-mode behavior throughout the bounce: see the main text. 
Figure (b) presents the time derivative of which clearly shows that the initially provided d^-mode is present even after 



the bounce, this coincides with the analytic solution in Eq. (|70[) . Meanwhile, Panel (a) shows that the initially introduced 
d a /3-mode is dominated by a newly generated C^-mode after the bounce. We take the initial condition in Table HVl thus k — 0. 
Even in the case of nonvanishing k, the large-scale condition which compares k 2 with a" /a term in Eq. (|64[) is well satisfied 
near the bounce. 
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FIG. 18: Panel (a) shows the behavior of initially near C Q/ 3-mode which transits to the d a) 3-mode after the bounce. As we are 
considering k = limit, the precise C^-mode in Fig. [16] sets C^l' = 0. In this Figure we introduced a small negative value of 



C aB= -1-9463109018142E-07X C^L This means that the present initial condition is a mixture of the dominating C a /3-mode 



-.(0 



plus a initially negligible but non-vanishing amount of the d a/ 3-mode. Panel (b) presents the time derivative of C^p which 
shows that the apparent d a 0-mode in the expanding phase is actually the one provided initially. 
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FIG. 19: Power spectra and typical evolutions for the C-mode of ip v . Panel (a) shows the initial and the final power spectra 
for our C-mode of ip v . Panel (a) shows that the spectral index is ns — 1 = 6 which agrees with our analytic result in Eq. 
(|76|) . Panel (c) shows the initial and the final power spectra for the analytic C-mode of (p v . In this case, only for large-scale 
(log k < —4.5) perturbation the spectral index coincides with our C-mode case. Panels (b) and (d) show typical evolutions of 
ip v corresponding to Panels (a) and (c), respectively, for two different wavenumbers. Panel (d) shows why the power spectrum 
at large wavenumber is enhanced compared with the one at small wavenumber. In Figs. I19l22l we take initial epoch at rj = —500 
[sign(a')-log(ai/a m i„) ~ —2.65] and final epoch at rj = 300 [sign(a')-log(a//a»n»n) ~ 2.25]. In order to get these power spectra, 
as the initial conditions both for our and analytic C- and d-modes, we use the Hankel function solutions in Eq. ([73} in the large 
scale limit. In order to find our C- and d-mode initial conditions for different wavenumbers we use the same method explained 
in Sec. lIVA3l 



[12] L.E. Allen, D. Wands, Phys. Rev. D 70, 063515 (2004). 
[13] F. Finelli, R. Brandenberger, Phys. Rev. D 65, 103522 
(2002). 

[14] D. Wands, Phys. Rev. D 60, 023507 (1999). 
[15] R.V. Buniy, S.D.H. Hsu, Phys. Lett. B 632, 543 (2006). 
[16] F. Lucchin, S. Matarrese, Phys. Rev. D 32, 1316 (1985). 
[17] J.M. Bardeen, Cosmology and Particle physics, edited by 



L Fang and A Zee (London:Gordon and Breach) (1988). 
[18] J. Hwang, Astrophys. J. 375 443 (1991). 
[19] J. Hwang, H. Noh, Phys. Rev. D 71, 063536 (2005). 
[20] J. Hwang, H. Noh, Class. Quant. Grav. 19, 527 (2002). 
[21] E.M. Lifshitz, J. Phys. (USSR) 10, 116 (1946); L.P. Gr- 

ishchuk, Sov. Phys.JETP 40, 409 (1975). 
[22] M. Abramowitz, I. A. Stegun (Eds.), Handbook of Math- 



34 



-20 



-40 



-60 



"i — i — i — | — i — i — i — | — i — i 1 1 r 

Our C-mode of <p 

initial power spectrum 

final power spectrum 



10 



-6 
log(k) 



cv -20 



-40 



-60 



I — i — i — | — i — i — i — | — i — i 1 1 r 

Analytic C-mode of ip 

initial power spectrum 

final power spectrum 



J I I I I I I I I L 



(c) 



10 -8 -6 -4 -2 

log(k) 



s. 
o 

o -5 

6 

i 

^ 1 n 
O - 10 



15 



1 1 1 1 1 1 

k = 


iii 

10- 3 


I I I I I I I I I I 


_ K = 


1 n-g 

1U D 


- 








: ^ 

1 1 1 1 1 1 1 i 1 1 


r- = 

fi ~ 

i i i i i i i i i 



(b) 



1 



sign(a,) -\og{cb/ch n 



10 



10 



t — i — I — i — i — i — i — I — i — i — i — i — r 

k 10 :i 



k=10- 6 




J I I I I L_ 



_L 



(d) 



-2 -1 
i 

sign(a,) -log(cfr/ct_ 



FIG. 20: The same as Fig. 1191 now for C-mode of ip x . Panel (a) shows that the spectral index for our C-mode is also ns — 1 = 6 
which agrees with our analytic result in Eq. (|76[) . In case of analytic C-mode, only for large-scale (logfc < —4) perturbations 
the spectral index coincides with our C-mode case. In this case, however, due to time dependent evolution of (p x , the amplitude 
of power spectrum simply depends on our choice of initial and final epochs. Panels (b) and (d) show typical evolution of ip v 
corresponding to Panels (a) and (c), respectively, for two different wavenumbers. Panel (d) shows why the power spectrum at 
large wavenumber is enhanced compared with the one at small wavenumber. 
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FIG. 21: The same as Fie. 1191 now for d-mode of ip v . Panel (a) shows that the spectral index for our d-mode gives near Harrison- 
Zel'dovich scale-invariant spectrum (ns — 1 = 0) which agrees with our analytic result in Eq. (I76|l . Panel (c) shows that the 
spectral index of the analytic d-mode also gives ns — 1 = 0. Panels (b) and (d) show typical evolution of (p v corresponding to 
Panels (a) and (c), respectively, for two different wavenumbers; these evolutions show why the final spectral indices resemble 
the initial ones. In both Panels (a) and (c), due to time dependent evolution of ip x , the amplitudes of power spectrum depend 
on our choice of initial and final epochs. For log(fc) > —2.5 the scale comes inside the horizon when the final spectra are 
measured. Although the figure is not presented in [13 |. based on the power spectrum in Panel (c) the authors of jl2j ] reported 
that the bouncing world model can generate scale-invariant spectrum. 
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FIG. 22: The same as Fig. 1191 now for d-mode of (p x . Panel (a) shows that the spectral index for our d-mode is ns — 1 = —4 
which agrees with our analytic result in Eq. (|76[) . In case of analytic d-mode, only for large-scale (log k < —4) perturbation 
the spectral index coincides with our d-mode case. In both Panels (a) and (c), due to time dependent evolution of ip x , the 
amplitudes of power spectrum depend on our choice of initial and final epochs. Panels (b) and (d) show typical evolution of 
ip v corresponding to Panels (a) and (c), respectively, for two different wavenumbers. Panel (d) shows why the power spectrum 
at large wavenumber is enhanced compared with the one at small wavenumber. 
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FIG. 23: This figure shows the final power spectra for ini- 
tial C a f3- and da/3- modes. The initial d Q( 3-mode has a scale- 
invariant spectrum with n T g — 0, and C^-mode has very 
blue spectrum with n Ti c a g = 6; these agree with our analytic 
result in Eq. (|80p . The amplitudes of power spectra depend 
on background properties. For log(fc) > —2 the scale comes 
inside the horizon when the final spectra are measured. 



